Solving the Cusp-Core Problem with a Novel Scalar Field Dark Matter 
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Matos, Guzman and Nunez (MGN) proposed a model of galactic halo based on an exponential- 
potential scalar field that could induce a rotation curve that is constant for all radii. We demonstrate 
that with suitable boundary conditions, such scalar field dark matter (SDM) can not only produce 
the observed constant rotation curve at large radii but also give rise to the correct power-law scaling 
near the galactic core region. This solves the existing cusp-core problem faced by the conventional 
cold dark matter (CDM) model. 
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I. INTRODUCTION 

The nature of dark matter has always been an im- 
portant issue in the study of cosmology. The Cold Dark 
Matter (CDM) model, which assumes the dark matter as 
weakly interacting massive particles (WIMPs), has been 
very successful in explaining the CMB anisotropies and 
the large scale structures in the Universe. However, it 
also has encountered several challenges such as the miss- 
ing satellites problem [IH2] and the cusp-core problem 

mm- 

The cusp-core problem arises from the mismatch be- 
tween the density profile obtained from the simulations 
based on the CDM model and the actual observations. In 
the core region of the galactic halos, the density profile 
has been found to follow the power law scaling with the 
galactocentric radius r, that is, p(r) oc r~ a |H|B]. The 
computer simulations based on the CDM model predict 
that a > 0.8 | EHU], while the observations suggest 
that a ~ 0.2 instead HH [TTj- (Note that some au- 
thors, based on CDM, found a ~ 0.7 instead [IB]-) In 
particular, the NFW profile [7J H] with a = 1 (cusp) at 
the center of the galactic halo has been found to agree 
with the result of CDM simulations [i HQ [H] . The 
pseudo-isothermal profile with a = (core) at the cen- 
ter, on the other hand, fits the observation data with 
a smaller v 2 than NFW does. Evidently, a core is pre- 
ferred over a cusp. This implies that if the dark matter 
is truly dust-like and if the CDM simulations are correct, 
then there would be an under-supply of gravity near the 
galactic center. 

One can entertain multiple explanations to this cusp- 
core problem. For example, at small scales the galac- 
tic dynamics may become so complex that the straight - 
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forward CDM ceases to be valid. That means it is likely 
that the current state of CDM simulations may not have 
included additional physics, such as baryon physics |19l - 
|2"5] , necessary for addressing this problem. On the other 
hand, the cusp-core problem may indeed be a signal that 
reveals the deficiency of CDM model, despite its great 
success at large scales. In this paper we choose to take 
the latter approach. 

Within the framework of the CDM model, the rotation 
curve is directly related to the density profile through the 
relation [To] . 

^Gp{r) = 2 V - d ^ + C-f- (1) 

One can see that under the CDM model, a power-law 
density profile implies a power-law rotation curve, i.e., 
v cx r@ , where a = 2/3 — 2. Therefore the mismatch 
of the density profiles in the CDM cusp-core problem is 
equivalent to the mismatch of the rotation curves, which 
is what one obtains directly from observations. Accord- 
ing to Eq.(l), there is a one-to-one correspondence be- 
tween the rotation curve and the density profile under 
CDM model. In principle, therefore, there exists a CDM 
density profile solution for every given rotation curve. 
However, computer simulations based on CDM model 
can only produce the cusp density profile, which does 
not conform with the observed rotation curve near the 
center. 

Several attempts have been made to build alternative 
DM models with non-zero pressure. Models consider- 
ing dark matter as an ideal fluid with isotropic pressure 
is known to be inconsistent with observations. Models 
with anisotropic pressure, therefore DM behaves as a 
nonideal fluid, have been investigated and shown to be 
self-consistent I26428j. It has been suggested that scalar 
fields can be a good candidate for this type of dark mat- 
ter. Several scalar field dark matter (SDM)models have 
been proposed [28lt33] , although their considerations are 
not same as ours. 
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The behavior of such scalar field dark matter, however, 
is greatly dependent on the form of its self-interaction 
potential. Among the various proposals, the exponential 
potential is often considered [33]. Matos, Guzman and 
Nunez (MGN), for example, showed that a consistent so- 
lution of Einstein equation giving rise to a flat rotation 
curve exists, if the self-interaction potential is of the form 
of an exponential function. Such solution reproduces the 
rotation curve in the outer region of the galactic halo. 
[28] 135] Matos et al. also showed that their model could 
well explain other observations such as that of type la su- 
pernovae [36] . The exponential potential scalar field may 
find its theoretical connections with, for example, the ra- 
dion or dilaton field in higher-dimensional gravity theo- 
ries [371 IMj • It can also be found in higher-order gravity 
theories [3HI3II]- In addition, the exponential-potential 
scalar field can be induced by non-perturbative effects 
such as gaugino condensation [41] . It is therefore quite 
natural to wonder, would such nonconventional SDM be 
able to solve the existing cusp-core problem associated 
with the conventional CDM model? 

In this paper, we consider a novel scalar field with a 
potential V(§) = — aexp(— as dark matter. This 
potential is essentially the same as that in MGN: V($>) — 
— aexp(— &$) except that we take the absolute value of 
the field in the exponent, so as to ensure the field is 
bounded from below in the system. We demonstrate that 
such a novel scalar field DM can readily reproduce the 
profile of the observed rotation curve from the galactic 
core to its outer region where the rotation curve is flat. 
Furthermore, the corresponding x 2 of our model is in 
general smaller than that of the NFW profile but com- 
parable to that of the pseudo-isothermal (ISO) profile. 
It is noteworthy that, by invoking only two parameters, 
which is the case in NFW and ISO, we are able to resolve 
the problem. That is, the exponential-potential SDM 
model is able to solve the existing cusp-core problem as- 
sociated with the conventional CDM model in a simple 
construction. In the CDM model the rotation curve can- 
not be uniquely determined owing to the degeneracy of 
the Einstein equation under zero pressure [5]. The de- 
termination of the rotation curve should therefore rely 
on simulations. In contrast, the rotation curve and the 
profile of energy-momentum tensor can be uniquely de- 
termined in our SDM model once the potential of the 
scalar field is given. Therefore, further simulation under 
our model is unnecessary once the unique solution that 
fits the observation is found. 

The organization of this paper is as follows. In Sec. II, 
we introduce the basic setting of the galactic model pro- 
posed by Matos, Guzman and Nunez. But instead of a 
constant, we invoke a hyperbolic tangent function for the 
rotation curve, which rather faithfully represents the ob- 
served generic rotation curve as a function of the galactic 
radius. This represents the first step of our attempt to 
solve the cusp-core problem in one stroke. We then show 
that the SDM with an exponential self-interaction poten- 
tial can indeed induce the desired rotation curve. In the 



Sec. Ill, we numerically compute the rotation curve in- 
duced by an exponential-potential scalar field dark mat- 
ter with suitable boundary conditions. The resulting 
rotation curves are shown to fit the observational data 
very well. We also include the baryonic energy density in 
the Einstein equation and show that such setting is self- 
consistent. In the last section we discuss the implications 
of our results. 



II. THE SELF-INTERACTION POTENTIAL 

Assuming spherical symmetry of galactic mass distri- 
bution, which includes both DM and baryons, we write 
down the general form of the metric: 



ds 2 



2 Mr)^2 + e 2X(r) dr 2 + ^2 + ^^2^2 



(2) 



We note that the actual mass distribution is not exactly 
spherically symmetric since the stellar distribution is of- 
ten non-spherical. Nevertheless, spherical symmetry has 
been shown to be a good approximation in the calcula- 
tion of the rotation curve [H] . To investigate the rotation 
curve governed by a scalar field DM, we write the Ein- 
stein equations as 



Ttot „ \ 



(3) 



where the total energy-momentum tress tensor is con- 
tributed by the SDM and the baryon matter: 



T 



T; 



SDM 



(4) 



With the aid of the Klein-Gordon equation that governs 
the SDM, 



dV 



0, (5) 
the energy-momentum tensor T^ M can be expressed as 

^,u-lg^ a ^,a-9„uV(^). (6) 



pSDM 



On the other hand, in our treatment of the baryonic 
energy- momentum tensor T* , which is predominantly 
contributed by the stellar mass, is an input that is deter- 
mined through the observation of the galactic luminosity. 
The Einstein equation then reads 

R^v = 87rG[$ + gn V V($) + baryonic terms]. (7) 

This is the only equation that governs our model. 

Under spherical symmetry, the g tt component of the 
metric is directly related to the rotation velocity of the 
galactic halo through the equation [26] 



ra(r)' 



c - 



(8) 
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From Eqs.(7) and (8), we see that once the rotation curve 
is determined the number of independent equations is the 
same as the number of undetermined variables: V($), 
<5(r), a(r) and A(r). 

In their model, Matos, Guzman and Nunez assumed 
the flatness of the rotation curve, with which they ob- 
tained the form of cr(r) through Eq.(8). Inserting it into 
the Einstein equation, they found V($) to be an expo- 
nential function |28j . 

To investigate the property of SDM near the galactic 
core, where the rotation curve is no longer flat, we replace 
the constant by a hyperbolic tangent function to model 
the rotation curve for the entire range of radius, i.e., 



v(r) — v n tanh(/3r a ), 



(9) 



where the radius r is in units of kpc. This rotation veloc- 
ity follows a power-law function at small r but saturates 
to a constant, vo, at large r. Such function fits the ro- 
tation curve reasonably well as shown in Fig. 1. When 
inserting the corresponding form of cr(r) in the Einstein 
equation, we keep only the leading order terms of <r(r), 
A(r) and v 2 /c 2 in order to streamline the calculation. We 
note that such simplification is reasonable since the ro- 
tation velocity is small compared with the speed of light 
and that the gravity inside the halo is weak [2H] ■ In par- 
ticular, it can be verified that the Matos-Guzman-Nunez 
solution remains unchanged under such approximation. 



where k — 1/8%G. These equations are still too tedious 
to render simple analytic solution. However if we solve 
them numerically, we see, as shown in Fig. 2, that the 
actual solution for the potential comes very close to an 
exponential function. This gives us confidence that the 
exponential potential may indeed generate the desired 
rotation curve under suitable boundary conditions. 




$ [10 24 eV] 



FIG. 2: Numerical solution for the potential V($) in Eq.(lO). 
The values of vo, d and a used in this calculation are taken 
from the fitting of the rotation curve data of NGC 6689. The 
dot-dashed curve is the numerical solution of V(<&) while the 
solid curve is the exponential fit. 
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FIG. 1: Using the hyperbolic tangent function to fit the dark 
matter contribution of the rotation curve of NGC 6689. 

Neglecting the baryonic contribution for the time be- 
ing, the Einstein equation reads: 

- - ^$^r a - 2 sech 2 (dr a )t & nh(dr a ) 
r c 2 

= ±[(l + 2A)V($) + i*'(r) 2 ] 

K I 

^ + ^ = I[(1 + 2A)F($) + V(r) 2 ] 

=^ + % tanh 2 (dr«) = + 2A)V(«) + ^'(rf), 

(10) 



III. CONSEQUENCES OF THE EXPONENTIAL 
POTENTIAL 



A. Rotation Curve 



Guided by the numerical solution, we now turn around 
and assume that the SDM potential is indeed exponen- 
tial: V = — aexp(— 6$), where a and b are positive. We 
substitute it into the Einstein equations and solve for the 
rotation curve. Ignoring the higher order terms as we did 
in Eq.(10), we have 

9 X 9 \' 1 1 

-^ + *-a" = V(*)(l-2A) + is'(r) a ], 
r r k I 

-^ + ^ = ^[-F($)(l-2A) + ^(r) 2 ], (11) 
where pb is the baryon density. 
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FIG. 3: Comparisons of the numerical solutions with the ob- 
servational data. The black dots are the observed rotation 
curves. The gray dots are the contribution of the dark mat- 
ter, and the dotted curve is the contribution of the stellar 
disk. The dashed curve is our result that includes only the 
dark matter, while the dot-dashed curve is our result with the 
baryonic contribution taken into account. For NGC 2976 the 
long dashed curve represents the contribution of HI and the 
solid curve represents the contribution of H2 



The task of obtaining an explicit solution to these 
equations is still formidable. We therefore resort again 
to the numerical solution. We observe that in every term 
of these equations the lowest order of the a derivatives 
is one. Thus we only need to specify the boundary con- 
ditions of A, a' and <&. We find that if the boundary 
conditions of all these three quantities are zero when r 
is extremely small, then the resulting rotation curve has 
the same shape as that obtained from observations. Fur- 
thermore, the resulting $(r) is always positive definite in 
our model, the same as in MGN. This indicates that our 
replacement of — aexp(— 6$) by — aexp(— does not 
alter any basic notions of MGN except to insure that the 
quantum states are bounded from below. 

Based on this solution, we analyze the five galaxies 
from Simon et al. [321 E3] an d two of the 19 galaxies 
from Blok et al. [45] whose Xnfw surpasses Xiso tne 
most, meaning that in these samples cores are greatly 
preferred over cusps. Among all of the analyzed galaxies, 
NGC 5963 is the only one with Xsdm larger than Xnfw 

Simon et al. [321 S3] have considered both the maxi- 
mum and the minimum disks. We note that in their anal- 
ysis of the maximum disk, they subtract the contribution 
of the stellar disk from the observed rotation curve before 
they fit it to the NFW and the pseudo-isothermal profiles. 
In their analysis of the minimum disk, however, they fit 
the observed rotation curves directly without subtract- 
ing the stellar disk contribution. Therefore we choose 
to compare our model with the other two mass models 
under the assumption of maximum disk only, which, as 
commented by Simon et al., is more physical. Before do- 
ing that we first use Eq.(8) and Eq.(ll) without including 
the baryonic density profile to calculate the SDM contri- 
bution alone. The comparison of Xsdm an d the x 2 of the 
other two profiles are shown in Table. 1, from which we 
see that except NGC 5963, all Xsdm are smaller than 
Xnfw an d comparable to Xiso- ^ n tne case °f NGC 
4605, Xsdm IS even smaller than Xiso- Comparisons of 
our numerical solutions to the rotation curve with the 
observational data are shown in Fig. 3. 

In our calculation of the total rotation 
curve, we include the stellar disk contribu- 
tion by invoking the fitting function, v^(r) — 

1.97 Wb 2 (r opt )(r/r opt ) 1 - 22 /[(^o P t) 2 + 0.78 2 ] 1 - 43 , which is 
commonly used [29, 46 . However the gas contribution 
differs from galaxy to galaxy, so in our actual calculation 
we used a four-to-six term polynomial to spline the 
discrete data points. We see from Fig. 3 that both 
results, the inclusion and the exclusion of the stellar 
disk contributions, fit the data reasonably well. Note 
that in order not to make the graphs too busy, only the 
error bars of the dark matter contribution are included. 
After all, that is what we need for the calculation of the 
Xsdm un( ier the maximum disk assumption. 

In the paper of Blok et al. [45], the data was ana- 
lyzed through a different method. They considered two 
cases where the 7*' 6 , that is, the light-to-mass ratio at 
3.6/^m, is either fixed through the mass model or kept 
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as a free parameter. They then considered two probable 
descriptions of the dark matter distribution, the NFW 
profile and the pseudo-isothermal profile. The total ro- 
tation curve, which is the square root of the quadratic 
sum of the rotation curves contributed by all the compo- 
nents, is then fitted to the observed rotation curve from 
which the corresponding \ 2 is calculated. 



TABLE I: x 2 of the data from Simon et al. under maximum 
disk. 
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FIG. 4: Comparisons of the numerical solutions with the ro- 
tation curves based on the observational data from Blok ct 
al. [45]. The black dots are the observed rotation curves, the 
dotted curve is the contribution of stellar disk, and the long 
dashed curve is the contribution of gas. The black curve is 
our result that takes the baryonic contribution into account. 
The dashed curve is our result that considers only the dark 
matter. 

For the two galaxies from Blok et al. [25], we assume 
that is the value that minimizes the Xiso an( ^ use 
Eq.(ll) and Eq.(8) to calculate the total rotation curve 
and the corresponding Xsdm- We find that although we 
do not treat 7^' 6 as a free parameter, we readily obtain 
a Xsdm that is much smaller than Xnfw- Comparisons 
of our numerical solutions to the rotation curves with 
the observational data from Blok et al. 45] are shown 
in Fig. 4. Once again, we find that our solutions fit the 
observations quite well. 



Energy-Momentum Tensor and Energy 
Conditions 



Having the desired rotation curve successfully deduced 
from the SDM exponential potential, we now inspect the 
property of the corresponding energy and pressure den- 
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TABLE II: x 2 of the data from Blok et al. 
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sity profiles of the SDM halo. The density and the pres- 
sure profiles of the dark matter can be expressed in terms 
of the potential and the field via Eq.(6): 

p d ar k = \e 2Mr) <5>'{r) 2 + V(r), 
Pr = \e 2X{r) <S>'{rf-V{r), 



Pi 



\e 2X(r) <S>'{r) 2 - V(r). (12) 



It is evident from these equations that the SDM model 
satisfies the non- ideal fluid equation of state, Pdark(f) = 
~Pt{r). 

We again take NGC 4605 as an example. Fig. 5 shows 
the energy density profiles, where pdark corresponds to 
the SDM halo and ptot includes both DM and the bary- 
onic matter. The radial and tangential components of 
the pressure are shown in Fig. 6. From these figures we 
see that the density and the pressures are of the same 
order of magnitude. This implies that the two terms, 
including the potential, on the RHS of Eq.(12) are also 
having the comparable magnitude. 

We see that although the energy density profile of the 
SDM is negative, the total energy density that includes 
the contribution from the baryons is always positive def- 
inite and therefore the system satisfies the weak energy 
condition [38] . We note that the dark matter sector it- 
self, though violating the weak energy condition, still sat- 
isfies the null and the strong energy conditions [55] ■ We 
therefore consider our solution reasonable. In addition to 
NGC 4605, we analyzed five additional objects: DD0154, 
NGC 2976, NGC 5949,NGC 5963, and NGC 6689. Again 
the null energy condition and the strong energy condition 
are satisfied in the region of interest, although the weak 
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FIG. 5: The solid is the density profile for the dark matter. 
The dotted curve is the density profile considering both the 
baryon and dark matter. 
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FIG. 6: The dotted curve is the profile of tangential pressure. 
The dashed curve is the profile corresponding to the radial 
part of the pressure. 



energy condition is not always satisfied. 

A negative density profile may seem unphysical at first 
sight. However, this is not true. In the conventional 
CDM model, owing to the fact that it provides no pres- 
sure, the energy density profile must be positive so as to 
render a positive, or attractive, gravity, which is the only 
required property for DM based on observations. For 
galactic models that are semi-Newtonian, this require- 
ment remains unchanged. However, in our model the 
magnitude of the pressure is comparable to that of the 
energy density. This changes the situation drastically. 
That is, in our case the strong energy condition ensures 
the attraction of the gravity without relying on the posi- 
tivity of the energy density, thanks to the nontrivial and 
substantial contribution of the pressure. 

To be more specific, it can be seen in the work of Faber 
and Visser [17] that the effective mass distributions as- 
sociated with the rotation curve and the gravitational 



m RC (r) = / drAirr 2 [p(r)+p r (r) + 2p t {r)] (13) 



and 



p{r) + -{ P r{r) + 2 Pt {r)) 



(14) 



respectively. In our model both [p(r) + p r (r) + 2p t (r)] 
and [p(r) + \{p r (r) + 2p t (r))] are positive definite for all 
r, which ensures that our effective mass distributions are 
positive and that the gravity is attractive. 

We note that negative energy density profiles do exist 
in physics. For example the energy density of Casmir 
energy is often negative. There also exist dark matter 
models that imply a negative energy density profile [551 
[29] . In the model of Arbey et al., a repulsive gravity in 
some region was even considered. 

In the conventional CDM model the rotation curve 
is solely governed by the DM density profile through 
Eq.(l). In our model where the pressure is nonzero and 
anisotropic, the energy-momentum that drives the stel- 
lar rotation velocity p(r) is replaced by Pdarkfr) +p r (r) + 
2p t (r) 47J. That is, 

^G[p dark (r) + p r (r) + 2p t (r)} = 2~ + (-) 2 . (15) 

r or r 

However, we also know that in our model Pdark{r) = 
—pt(r). This means that the rotation curve is determined 
by the sum of the two pressure components, p r +pt (but 
not the total pressure P tot 



■pI 



however). It is 
interesting to note that while in CDM the rotation curve 
is dictated entirely by the density, in contrast in SDM it 
is entirely governed by the pressure. As for the relative 
contributions between the two pressure components, we 
see from Fig. 5 that both the radial and the tangential 
pressures exhibit similar radial dependence with roughly 
the same magnitude, except that the radial component is 
consistently a few times larger than that of the tangential 
one over essentially the entire range of the radius. 

Evidently the derivative of our potential, 
— aexp(— b\Q\), is discontinuous at $ = and the 
model is therefore ill-defined. On the other hand, we 
found that the resulting rotation curve in our model is 
totally insensitive to the the profile of potential around 
$ = 0. One may therefore regard our potential as a 
phenomenological model that is valid for all values of 
$ except 0. We assume that the true potential is the 
combination of this exponential potential and a small 
correction term that rounds off its sharp edge at <f> = 0. 
Being minute and effective only at $ ~ 0, such addition 
to the potential would not affect any of our profile 
analysis. 



IV. CONCLUSION 

We have shown that the scalar field dark matter with 
an exponential self-interacting potential can successfully 
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generate a unified and universal stellar rotation curve 
that covers from the galactic core to its outer region. 
Moreover, the Xsdm based on our scalar dark matter 
model is found to be generally comparable with xjso 
and smaller than Xnfw- f n sorne of the galaxies, we 
even found Xsdm smaller than xjso- The only galaxy 
with Xsdm > Xnfw i- e -' NGC 5963, may be an excep- 
tional case as commented by Blok et al. The brightness 
of its inner disk indicates that it may not be dark mat- 
ter dominated all the way to the center [6]. We recall 
that the pseudo-isothermal profile (ISO) agrees with the 
observation data whereas the NFW profile, which con- 
forms with the CDM-based simulations, does not agree 
with the observation data near the galactic center. Being 
more consistent with the ISO profile and thus the obser- 
vations, the scalar dark matter (SDM) model ameliorates 
the cusp-core problem associated with the CDM model. 

In the conventional CDM model, there are in general 
three unknowns to be solved in the Einstein equations un- 
der the assumption of spherical symmetry. These are the 
DM density profile Pdark(r) and the two metrics gtt(r) 
and g rr (r). However, since the CDM has no pressure, 
the Einstein equations are degenerate and that renders 
only two independent equations in the leading order ex- 
pansion |26j . As a consequence, only the relationship 
between these three unknowns, instead of their absolute 
values, can be deduced. This is why it is not possible to 
derive the rotation curve from CDM model based on first 
principles. One usually relies on numerical simulations 
to discriminate possible solutions [B] . In our SDM model 
there are also three unknown quantities, <E>(r), g tt (r) and 
g rr {r), to be determined after the form of V($) is speci- 
fied. However the difference is that in our case the pres- 
sure is nonzero and thus all three Einstein equations are 
non-degenerate. The three unknowns can then be deter- 
mined explicitly The profiles based on the SDM model 
are thus attainable from first principles without the need 
to resort to simulations. The core (instead of cusp) pro- 
file so deduced from SDM is therefore unambiguous. We 
find this aspect rather novel. 

We have demonstrated the universality of the form of 
our SDM potential among the galaxies that we investi- 
gated, with a ~ 0(l(r 2 )eV and b ~ 0(l(r 24 )(eV)-\ al- 
though these values differ slightly from galaxy to galaxy. 



Comparing these values with that in MGN [2S], we find 
that the exponent of our potential comes very close to 
that in their model, where b ~ 2\J kc 2 Jvq. Such cor- 
respondence is expected. MGN considered only the case 
where the rotation curve is a constant, which corresponds 
to our model at large r. If the rotation curve is constant 
for all r, the Einstein equation is analytically solvable, 
which has been carried out in MGN's work. Since our 
result should approach that of MGN at large r, the value 
of our b is expected to be similar to theirs. 

Although we mentioned in the introduction that the 
exponential potential is quite common in higher dimen- 
sional or higher order gravity theories, so far we are not 
able to deduce the coefficients in the exponential poten- 
tial from first principles. Our model should therefore be 
viewed as a phenomcnological or effective theory for dark 
matter halo. Even so, from the phenomenological point 
of view it is remarkable that by replying on no more pa- 
rameters than that in NFW, our model can solve the 
cusp-core problem. 

It is encouraging that our novel scalar field dark matter 
model with an exponential potential is able to resolve the 
cusp-core problem that stems from the cold dark matter 
model. Before this SDM model can be taken more seri- 
ously, however, it is essential that this non-conventional 
SDM can reproduce the large scale structure formation 
history of the universe, which is what the conventional 
SDM does its best. Another challenge is the need of an 
explicit SDM production scenario in the early universe 
based on a fundamental theory that would give rise to 
the desired amount of such SDM that would saturate the 
DM content of the universe. 
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